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Abstract. The asymptotic expansion method is generalized from the periodic set- 
ting to stationary ergodic stochastic geometries. This will demonstrate that results 
from periodic asymptotic expansion also apply to non-periodic structures of a certain 
class. In particular, the article adresses non-mathematicians who are familiar with 
asymptotic expansion and aims at introducing them to stochastic homogenization in 
a simple way. The basic ideas of the generalization can be formulated in simple terms, 
which is basically due to recent advances in mathematical stochastic homogenization. 
After a short and formal introduction of stochastic geometry, calculations in the sto- 
chastic case will be formulated in a way that they will not look different from the 
periodic setting. To demonstrate that, the method will be applied to diffusion with 
and without microscopic nonlinear boundary conditions and to porous media flow. 
Some examples of stochastic geometries will be given. 



1. Introduction 

Homogenization has become an important modeling tool for multiscale problems, i.e. 
phenomena that have causes and effects on multiple spatial scales. The applica- 
tion ranges from Physics (porous media flow, freezing processes in porous media, 
erosion[2, 7, 22j), engineering (composite materials, reaction diffusion equations in 
catalysts [19t .llj ) to biology (processes in tissue||13j, membranes[20j). Also the numer- 
ical investigation of such problems is of interest (see jl6] and references therein). 

Homogenization considers physical, chemical or biological processes on large domains 
with a periodic microstructure of period £ <^ 1. On this periodic structure, a set 
of equations describing the several processes are set up and the limit behavior of the 
solutions of these equations is investigated as e — )■ 0. Note that e — )■ is a mathematical 
equivalent of the assumption e ^ 1, which means basically very small microstructures 
compared to the domain of interest. For example, one may think of a solid structure 
perforated by a periodic structure of channels, which are filled by a Newtonian fluid, i.e. 
a fluid whose motion is described by the Navier-Stokes equations. Since one is usually 
interested in a scale of meters and the pores have sizes of mm, it follows e ^ 10~^. In 



the limit e — )■ 0, the velocity field would follow Darcy's law [lOj (see sections 2.2 and 



4.2). 



There are two ways to obtain the limit equations: Via strong mathematical calculations 
(proofs of convergence of solutions) or via rather formal calculations. One of the formal 
ways to obtain averaged equations is the so called asymptotic expansion. Its ansatz 
is to expand the solution in a series of functions multiplied with increasing powers of 
e, where the functions in this expansion depend on a global variable and a periodic 
variable. We will demonstrate the basics of this approach below in section [2j 

Strong mathematical investigations where originally based on strong and weak conver- 
gence methods in Sobolev spaces, which is nowadays still a method of choice. However, 
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lots of other convergence methods entered homogenization theory such as T- and G- 
convergence, see |26j for an overview. 

In the beginning of the 1990's, Allaire and Nguetseng [21], [1] developed the so called 
two-scale convergence which was later extended by Neuss-Radu, Zhikov and Lukkassen 
and Wall [T8| [25l [T2] . As an interesting feature, two-scale convergence is related to 
an asymptotic expansion in the weak sense and therefore closely connected to the 
formal method. This similarity becomes even more striking by the method of periodic 
unfolding, developed by Cioranescu, Damlamian and Griso [5i 14], and which can be 
considered as a consequent generalization of two-scale convergence. In comparing the 
methods, one may get the idea that every result from formal asymptotic expansion 
could be proved rigorously. However, for many homogenization problems there are lots 
of technical difficulties like regularity proofs or a lack of Poincare-inequalities which 
are sometimes to hard to be overcome. 

There is, however, the justified criticism, that the methods above are based on the 
assumption of a periodic structure and that the resulting equations therefore may 
not be valid for a non-periodic medium. Nevertheless, we expect that many of the 
derived equations may also hold for non-periodic microscopic geometries: Darcy's law is 
observed to hold in most porous media may it be sand, silt or loam. In a polycrystal or a 
composite material (industrial ceramics) with their complex microstructure, we assume 
that heat transport still follows Fourier's law with averaged coefficients. Reaction 
diffusion equations derived for a catalyzer should also hold if the microscopic structure 
is not perfectly periodic. 

So, it is an important question whether or not it is possible to apply homogenized 
results to non-periodic geometries and under what circumstances. Mathematicians 
were aware of this problem and tried to apply homogenization techniques to non- 
periodic structures. An overview over many results before 1994 can be found in [26j. 

In 2006, Zhikov and Piatnitsky [27] where able to introduce a two-scale convergence 
method on stationary, ergodic stochastic geometries on a compact probability space. 
Their results where generalized to arbitrary probability spaces by the author in [8j . A 
former attempt by Bourgeat, Mikelic and Wright [3j in the non-ergodic setting used 
an averaged form of two-scale convergence and was not applicable to problems on 
manifolds or complex boundary conditions. Also, due to the averaging, it was less 
close to the periodic setting than [27] . 

The results in [8] provide a geometrical interpretation of subsets of the probability 
space and demonstrates that all periodic quantities find their precise analogue in the 
stochastic case and vice versa. Moreover, if a stationary ergodic structure is periodic, 
the probability space has to be essentially the unit cell. The theory has been applied 
successfully to heat transfer in polycrystals in |9j. 

Due to the connection between two-scale convergence and asymptotic expansion in 
the periodic setting, these mathematical results give the basic idea to extend the as- 
ymptotic expansion to the stochastic case. Note that it is not the intention of this 
article, to give a rigorous mathematical introduction to stochastic geometry or sto- 
chastic homogenization. Rather it is the aim of this paper to demonstrate that it can 
be mathematically justified to apply results from asymptotic expansion to stochastic 
settings and that there is formally no difference in the calculations. The most important 
theoretical results are cited in the appendix. For proofs of the fundamental theorems 
which are cited in this article, the reader is referred to |27l i8j and the references therein. 

This article is organized as follows: In the next section, the basic idea of periodic 
asymptotic expansion is explained and the resulting equations for two sample problems 
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(diffusion and porous media ffow) are given. In section |3| the concepts of stochastic 
geometry will be introduced in a formal and mathematical non rigorous way. The 
connections to periodic geometries are explained. In section|4]the asymptotic expansion 
on stochastic geometries will be introduced and applied to thermal diffusion, porous 
media flow and diffusion processes with reactions on the microscopic boundary. In 
section [5] some simple examples of stochastic geometries will be given. 



2. Recapitulation of the Periodic Case 

In this section we will shortly recapitulate the asymptotic expansion technique for 
periodic structures. It is not the intention of this section to go into the details of 
calculations but to rather explain what are the mathematical problems, the ansatz 
that people use to solve these problems and what the results look like. For both 
examples, precise calculations will follow in section |4] below. We will first consider the 
standard homogenization problem of diffusion (following [lOj) and then discuss Navier- 
Stokes flow through porous media. In all calculations below and throughout the paper, 
n denotes the dimension of space (in a physical setting = 2 or n = 3). 

2.1. The standard problem. Given an open and bounded set Q C investigate 
diffusion processes (for example heat diffusion) with a source term and rapidly oscillat- 
ing coefficients which are due to a periodic microstructure. In particular, assume that 
for Y := [0, 1[" we have Y = AUBUT with A, B open in Y such that A fl 5 = and 
r = OA = dB in Y. We expand Y periodically to M" and identify A, B and F with 
their periodic continuation. We furthermore denote A^ = eA, B^ := eB and F^ := eT. 
The diffusion constants in A and B are thought to be different. 

The mathematical problem reads 

dtu' - div (L)^Vm^) = / on ]0, T] X g 

M^ = on]0,r]xaQ 
M^(0, ■) = a(-) on Q for t = 

where is given by := D{t,x, |) and D is a 1^-periodic symmetric tensor with 
bounded entries, while / is a function on Q which may also depend on time. As 
an example, we could have = 1 on A^ and = 10 on the interior of 5^ with 
D{t,x,y) = 1 + 9xb(z/). 

For small enough microstructures, the oscillation in D may have a minor effect on 
and the system may be described by some averaged diffusion coefficient _D^°™. It 
may therefore be sufficient to solve an approximating system with this smoothed 1)^°"^ 
without the strong oscillations of Z)^. Therefore, we investigate as e — by an 
ansatz 

(2.1) u^{t,x) = ^e'ui{t,x,-) 

i>0 ^ 

which is called the asymptotic expansion of and 

: ]R>o X Q X Y ^ M 

which are 1^-periodic functions in the third coordinate. The gradient operator turns 
into V = + -Vy while the divergence reads div = diva^ + -div^. Altogether, the 
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expansions above inserted into the diffusion equation yield 

dtuo - ^divy {D^VyUo) - Jdiv^ (D^VyUo) - Jdivj, {D^'VxUq) - Jdiv^ (D^VyUi) 

- div, (D'V^Uo) - div, (D'VyU^) - divj, {D^'^ ^u^) - div, {D'^ yU-,) + 0{e) = f 

Assuming tliat the terms of order e"^ and vanish, it will be shown below in section 
4.1 that this system simplifies to a single equation for mq 

with defined in an appropriate way. 

2.2. Porous media flow. Suppose that A"^ is connected and consider the Navier- 
Stokes equations on Q (1 A"^: 

dtv' + {v' ■V)v' ~ div {uVv') + Vp' = f on ]0, T] x (Q n A') 

divi;" = on ]0,T] X (Qn A") 

v' = on]0,T] X 

v' = on]0,T]xB' 

v%0, ■) = a{-, -) on Q for t = 

where / is an external force and a is 1^ periodic in the second argument. For an ansatz 

{x,y) ^ v^{x,y) {x,y) pi{x,y) 

such that the solution and can be described by 



oo oo 

X. ^ sr-^ 4 , X, 



^ e%i{x, -) f = ^ e'p-iix, -) 

the resulting set of equations is 

- dwy{uVyV2) + VxPo + VyPi = g 
(2.2) divyi;2 = 

V2{x, ■) = on dVi 



We will see in section |4.2| that for a suitable choice of a matrix K, V2 and po fulfill the 
following equation: 

/ V2 = K{g- Vr^po) 

J V 



which is Darcy's law. 



3. Stochastic Geometry 



In order to develop the asymptotic expansion method for the stochastic (non-periodic) 
case, it is necessary to develop a suitable mathematical formalism. This will be the aim 
of the first part of this section. As a necessary condition, we expect that the periodic 
case is but a specialization of the stochastic one, which will be shown in the second 
part. The basic idea of the theory described below is to replace 1^ by a probability 
space Q. It is then necessary to identify equivalents of A, B and F on similar to the 
periodic case. After that, the relations to periodic structures are pointed out, before 
going on the asymptotic expansion and examples for stochastic geometries. 
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3.1. A phenomenological introduction to stochastic geometry. Instead of the 
periodic cell Y, consider a probability space {Q,a,fi) with the probability set Q, the 
sigma algebra a and the probability measure fi. Assume that i7 is a metric spac^ and 
there is a family {tx)x&R" of measurable bijective mappings r^. : f2 h-)- which satisfy 



• A : M" X r2 — > i7 (x, w) (— 7- TxOJ is continuous 

The family is then called a Dynamical System. Additionally, we claim ergodicity of 
r,, which is one of the following two equivalent conditions 

[/(w) = fir^uj) Va; G M" , a.e. u eVL]^ [fi!^) = const for fi — a.e. u E Q] 
^^■^^ [P {{TxiB) U B)\irxiB) n B)) = 0] ^ [P(fi) G {0, 1}] . 

This condition seems to be only technical, but in fact, it is crucial for mathematical 
homogenization in [HI [27] and replaces periodicity as well as the requirement that Y is 
simply connected. 

In order to introduce stochastic geometries from a phenomenological point of view, 
assume that there is a measurable set A C such that for the characteristic function 
Xa '■ ^ ^ {0, 1} holds xa{t»^) '■ IR" — ^ {0, 1} is the characteristic function of a closed 
set A{u)) C M" yU— almost surely in u. A{u) is then called a random closed set|^ Note 
that it possesses the property Xa{uj) + y) = Xa^t^uj) in) which is called stationarity. 

For some of the examples below, we assume at the same time that there is B G Q such 
that B{u) = ]R"\y4(a;) is the closure of the complement of A{u). The set T := A (1 B 
is then also measurable and has the property that r(ci;) = dA{u). Indeed, as shown in 
the appendix, F can be considered as an abstract manifold, since it can be assigned a 
Hausdorff measure and a normal field. 

Since the space Q is metric, there is a set of continuous and bounded functions Cfy{Q). 
At the same time, since A and B are measurable subsets, it is also possible to consider 
continuous functions Cb{A) and Cb{B) on A and B. With help of the dynamical system 
r,, it is possible to also define derivatives according to 

(3.2) du,,if{uJ) ■= lim 



>o h 

where (ei)i<i<„ is the canonical Basis on M"and / G Cfe(fi). The space of continuously 
differentiable functions is then 



C\n) := {/ G C{n) I limit Q is defined Vw G n and d^^if G \/l < i < n} . 



Similar spaces may also be defined on A and B. We may also define the gradient V^; 
and the divergence div^ accordingly. For any set X, any function / : f2 — )■ X and 
any u E Q, the function f{T,uj) : M" — X is called a realization (w-realization) of 
/ and for any continuous or differentiable function the realization is also continuous 
or differentiable, respectively. This is due to the continuity of r. In particular, for 
/ G Cb{A) holds /(r.w) G Ci,{A{u:)) and similarly for B and F. 



""^which means there is a distance function d : x — !• R>o such that d(aji,a;2) = oji = ^2, 
d{uii,ui2) = d{ui2,LiJi) and d{u:i,uiz) < d{uii,LU2) + d{ijj2,ui^) for any uii,ui2,L03 G fi. But this is only 
technical to be able to use the notion continuity. 

^Mathematically, as shown in the appendix, a random closed set is defined via a mapping A{lo) into 
the set of closed sets. The existence of sets A,B,T such that the above properties are fulfilled is 
shown afterward. For formal calculations, it seems more appropriate to start the other way round. 
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Finally, for a scaled random set A'^{u) := eA{uj) holds 

X 

XAHi^)i^) = XA{w){-) = XA(l-fW) , 

and for a function / G Cb{A) holds f^nuj) G Cb{A'^{uj)) for all the realizations u. 

3.2. The periodic case. It is helpful to compare the results above to the periodic 
case. First, note that Y = [0, together with the Borel sigma algebra and the 
standard Lebesgue measure C can be considered as a probability space {Y, as, C). For 
any x G MJ^ let [x] G denote the vector of integers such that x G [x] + 1^. Define 
the following family of mappings 

T,:Y^Y 

y\-^y + x— [y + x] 

which has all the properties we claimed for a dynamical system to hold. For any closed 
set A gY, A{y) is automatically closed for all y G 1^ and A{0) is the periodic contin- 
uation of A. The continuous functions on Y coincide with the 1^-periodic continuous 
functions on M" and the derivatives dy^i coincide with the classical derivatives di. B is 
the closed complement of A in 1^ and F is the boundary of A in Y. 

It is therefore clear, that any method which is based on the formalism introduced above 
is at least applicable to the periodic setting. Section |5] will demonstrate that a much 
broader class of geometries is covered by this approach. 



4. Asymptotic expansion on stochastic geometries 

We are now able to introduce asymptotic expansion on stochastic geometries. The 
method will be introduced via a sample calculation for the standard homogenization 
problem, which is diffusion with rapidly oscillating coefficients. This will be done using 



an expansion (4.1) of the unknown similar to (2.1) in section 2.1 Afterward, the 
method will be applied to porous media flow and diffusion with nonlinear microscopic 
boundary conditions. The formal calculations are well known for the periodic case and 
most of them can be found in [lOj among others. Here, they will be presented in the 
stochastic framework to demonstrate the similarity with the periodic case. 

4.1. The standard homogenization problem. The standard problem in the sto- 
chastic setting reads 

dtu' - div {D'^Vu') = f on ]0, T] X Q 

u' = on]0,T]xdQ 
u%0, ■) = a(-) on Q for t = 

where now = D{t,x,T^uj). An example would be D{-,-,u) = 1 for a; G A and 
D{-,-,uj) = 2ioT u e B. 



While the original idea is to expand the functions m*^ asymptotically by (2.1) with 
functions Ui : ]R>o x Q x 1^ — )■ M, it will now be based on a stochastic ansatz. In 
particular, for a given choice of the microscopic geometry u with the micro structures 
^4*^(0;), B^iuj) and r^(a;), we may expand by 

(4.1) x) = e'^Ujit, x, t?luj) . 



ASYMPTOTIC EXPANSION ON STOCHASTIC GEOMETRIES 



The gradient and the divergence have the expansions 

(4.2) V = + ^V^, div = div^ + ^div^ 



Using an expansion (4.1) for with (4.2) will lead to 

dtuo - \diy^ {D'^V^uq) - -div^. {D^V^uq) - -div^ {D^\J ^u^) - -div^ (D^V^Mi) 

E'^ e e e 

- div,. (D^V.Mo) - div,. (Z^^V^Mi) - div^ (D^V,Mi) - div^ (D"V^M2) + 0{e) = f 
which means we split up the equation in terms 

£-2(...)+£-l(...)+£°(...) + 0(£) =0. 

Since the latter equation should hold for all e, it follows 

div^ {D{t, X, uj)'VujUo) = . 

Note that in the periodic case, the latter equation implies uo{y) = const, while in the 
stochastic case, this is not clear. However, due to Gauss' theorem [6] for dynamical 
systems it follows by testing the equation with uq'- 



/ D{t, X, u) I V(^'Uo(t, x,u))f dfi{u) = 
Jn 



which is VujUq = 0. Definition (3.1 )i of ergodicity yields uo{u}) = const. The terms of 
order yield 

div^ {D'V^Uo) = -div^ {D'V^u,) . 
Now, if is a solution to the cell problem 

div^ {D{t, X, uj)Vuj4>j) = -divuj {D{t, x, uj)ej) 
for j = 1, . . . , n, the function Ui can be expressed by 



a; 



Ui{t,x,uj) = (pjjt, X, uj)djUo{t, 

Existence of (pj with (pj = can be shown with help of the Poincare inequalities 
from the appendix and the Lax-Milgram theorem. 

Finally, the zero order terms add up to 

(4.3) dtUo - div,^ {DV^Ui - DV ^U2) - div, {DV ^Uq) - div, {DV^Ui) = f . 

By integrating the latter equation over Q with help of Gauss' theorem [6| one obtains 



9iMo-div, / L'rf/iV,Mo - / div, (DV^ni)= / /rf/x , 
\Jn J Jo, 

The third term on the lefthand side can be reformulated to 



which finally yields 



( div, {DVu>ui) = div, (^j '^uj(pjdjUo^ 

5iMo - div, (D"°"^V,Mo) = //rf/i. 

Jn 



Here, 



Jn 

is a symmetric positive definite matrix. The proof is the same as in [10 
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Remark. The obtained limit problem is independent on the particular choice of the 
realization oj which we used for homogenization. This is a feature of the ergodicity 
and stationarity and reflects our expectation that the averaged equations should not 
depend on the specific geometry but on "the type" of geometry. We will see that the 
other examples share this property of the limit equations. 

4.2. Porous media flow. Suppose that A"^ is connected and consider the Navier- 
Stokes equations on Q (1 A^: 





= f 


on ]0,T] X (QnA^u)) 




= 


on ]0,T] X (QnA'iuj)) 




= 


on ]0,r] X {dQur{uj)) 


v' 


= 


on ]0,r] X B'{io) 




= «(■,-) 


on Q for t = 



where / is an external force and a is 1^ periodic in the second argument. From physical 
investigations, we know that the flow through a porous medium obeys Darcy's law and 
it is the aim of the following calculations to demonstrate that it is possible to obtain 
it as a limit problem for vanishing porescale in the stochastic setting. 

For an ansatz 



Vi 



Pi : M>o xQ xn^R^ 

(t, X, u) I— i- pi(t, X, co) 

such that the solution and p^ can be described by 



(t, X, u) I— !■ Vi(t, X, co) 



(4.4) 



e^Vi{t, X, r|w) 

i=0 



00 



1=0 



the resulting set of expanded equations is up to order 0: 



(4.5a) 



(4.5b) 
(4.5c) 



'div^{uVu^Vo)) + e"^ {{vq ■ V^) Vq - div^{uV^Vi) - 2div^(z/V^t;o) + V^po) 
+ £° ((do ■ Vo + (vi ■ V^) Vo + {vo ■ V^) Vi) 
+ eo (-div^(z/V^t;2) - 2div^.(z/V^Do) + V^po + V^pi - g) + 0{e) = 

e~^div^vo + div^i^o + div^i;i + e (diva;t;i + div^i;2) + 0{e'^) = 

^£:*t)j(x, ■) = on dYi, ^ £:*i;i(-, w) = on dil 



For each power of e, a set of equations is obtained, which has to hold independently on 
all the other equations such that the whole group of equations is valid for all choices 
of e. 



The order of —2 in (4.5a) together with the order —1 in (4.5b), and (4.5c) yields for vq 

: —div^{i/'VuiVo) = on A 
(4.6) div^wo = on A 

Vo{x, ■) = on r . 
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Again, in the periodic case, it would immediately follow Vq = 0. However, Gauss' 
theorem |6] again provides the necessary framework since 

2 



implies WuiVq = which is together with ergodicity (3.1) Vo{u) = const and with (4.6)3 
finally vo = 0. 

The latter result together with the order —1 in (4.5a), order in (4.5b) and (4.5c) 
yields for vi 

- div^(z/V^i^i) + V^Po = 
(4.7) div^vi = 

vi{x, ■) = on dYi 



which is again vi = due to vi ■ V^jPo = and (4.7)2. 

Using these results in the zero-order term in (4.5a), order 1 in (4.5b) and order 2 in 
(4.5c), we get 

-div^{uVu;V2) + V^po + V^pi = g 
(4.8) div^t;2 = 

V2{x, ■) = on r 
Assuming that there are solutions to the problems 

-div^(z/Va,Mi) + V^^IIj = Ci 
div u,Ui = 
Ui{x, ■) = on r 

where Cj is the i-th coordinate vector of M.^, it is easy to see that there is pi and 
V2 '-=^{9 — ^xPQ)iUi such that {v2,Po,Pi) is a solution to (4.8). Existence of Ui may 
be again shown similarly to the previous example. Defining a matrix K by 



one may check that we obtain indeed Darcy's law which was expected: 

V2 = K{g- V^.po) 



L 



4.3. Diffusion with nonlinear boundary conditions. Based on the calculations 
in section 4J^ consider the following problem 

dtu' - div {DlVu') = f 
= 

M^(0, ■) = a(-, T'Uj) 
combined with the nonlinear boundary problem 

dtU' + g{u' = Q 

U'{0,.) = U{.,T,oj) 

This system describes diffusion with a production or uptake due to reactions on the 
microscopic boundaries. The e-factor in the boundary condition takes into account 
that the microscopic surfaces increase with a factor 6~^. From physical perspective, we 



on ]0,T] X {QnA%u)) 
on ]0,r] X dQ 

on]o,r] X (r(w)nQ) 

onQnA'{u) for t = 

on]0,T] X {r{uj)nQ) 
on P(a;)nQfor t = 0. 
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expect that the reactions on the walls will lead to a macroscopic production or uptake 
term in the limit problem. 



Using an expansion (4.1 ) for u'^ and a similar expansion for f/*^, one would immediately 
conclude for wq 

dtWo + g{uM) = ^ on]0,r]xQxr 

f/o(0,.) = ongxrfort = 0. 

For the expansion of one again obtains 

dtu^ - ^div^ {D'V^u^) - -div, {D'V^u^) - -div^ {D'V ,u^) - -div^ {D^V^u^) 
- div, (/^^V,.Mo) - div, (Z^'V^Mi) - div^ (D"V,Mi) - div^ (D'V^Ms) + 0{e) = f 
together with the following boundary conditions on F: 

(4.9) \ (DV^Uo) ■ur + -D (V^Uq + V^Mi) -ur + D ( V,Mi + V^u^) ■ z^r -^K, Uo) = . 
E'^ e 

It is again possible to obtain with slightly modified argumentation in the partial inte- 
grations Uq{uj) = const as well as 

n 

ui{t, X, u) = (pjit, X, uj)djUo{t, x) , 
i=i 

where now the (pj are solutions to 

div^ {D{t, X, uj)Vu}(pj) = — div^ x, oj)ej) on A 

{DV^(f)j + Dej)-ur = on F . 



The zero order terms in (4.3) are now integrated over A which yields together with the 



zero order boundary condition in ( |4.9[ ): 

5iMo - div4D'^°"^V.Mo) - / g{uo,Uo)di^= f fd^. 

Jr Jn 



where now 



J A 



We therefore find the nice property that the macroscopic unknown Uo is produced or 
consumed by microscopic reactions on the surfaces, which was expected. 



5. Non-periodic models captured by stochastic geometries 

We will now describe some basic modeling tools for stochastic geometries and give some 
simple but interesting models for different applications in microstructures. Note that 
for many real world applications, it is up to now not possible to give a model for the 
natural geometries. For example natural soil with roots, wormholes and other obstacles 
is currently out of reach. Therefore, the models given here can only be considered as 
very rough models for porous media and other applications. Also, we will restrict on 
the most simple models. For more models in stochastic geometry, refer to Stoyan, 
Kendall and Mecke |23j. 
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Figure 5.1. Left: Poisson-Voronoi tessellation with its underlying Pois- 
son point process; Right: The corresponding Delaunay tessellation. The 
point processes and tessellations where constructed using the "Stochastic 
Geometry 4.1" software, developed at the TU Bergakademie Freiberg, 
Institut fiir Stochastik. 



5.1. Point processes. A point process (PP) is a measurable mapping A : Q ^ (M") . 
Let V{N,Q,p) denote the probability to find points in the bounded and open set 
Q C M" for the point process p. For some open set C with V 3 0, define the 
intensity Ap by 

Ap:= hm -^EiV{;tV,p)) , 

where E{V{. . .)) is the expectation value of A^. According to |2S1 a PP is stationary 
if its characteristic is invariant under translation, i.e. V{N, Q,p) = V{N, Q + x,p) for 
all X G M". From a stochastic point processes, one may construct higher dimensional 
structures in a deterministic way. Such structures would still be stationary random 
closed set. We will come to that point below. 

A prominent example of a PP is the so called stationary Poisson point process which 
has the property that 

V{N,Q,p) = ^^^exp{-X,\Q\) 

An other useful class of PP are so called hard core PP: From a given random PP any 
point is erased if its nearest neighbor has a distance less than a certain value d. 

5.2. Voronoi- Tessellations and Delaunay-Diagram. Starting from a PP, one may 
construct the so called Voronoi-tessellation as the set of all points in M", whose nearest 
two neighbors of the PP have the same distance. This can be used for example as a 
model of polycrystalslSj: The point process represents the set of crystallization nuclei 
who start to grow at the same time with the same speed and who stop growing the 
moment they hit each other. 

A generalization of this model which models crystal growth starting at different times 
with different speed is the so called Johnson-Mehl-tessellation. 

As a dual of the Voronoi-tessellation, one can consider the Delaunay-diagram which 
connects all points of the PP who share a common border in the Voronoi-tessellation. 
Connecting the points with cylinders (pipes) instead of lines, one would already get a 
model for a porous medium. Examples of both are shown in figure |5.1| 
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Figure 5.2. A ball process based on a Mattern hard core point process 
of type 2 (generated with Stochastic Geometry 4.1) 



5.3. Ball- and Grain-models. Based on a hard core PP, assign to each point of the 
PP a ball with a fixed radius (Figure 5.2). This model could be used to model a porous 
medium in 2D. However, in 3D we expect the grains of the porous medium matrix to 
touch some neighbors and that none of them is free without touching any neighbor. 
To achieve such a geometry, consider the following construction: 

From a Voronoi-tessellation consider for each point Pc of the PP the set of points Pf, i 
where the corresponding Delaunay-diagram hits the grain boundaries. We may then 
interpolate these points by a smooth manifold to mark the surface of a grain. As an 
example, define 

The sets G{-) could then model the grains of sand and the complement of the G would 
be the porous medium. 



6. Conclusion 



We saw that asymptotic expansion can also be applied to stochastic geometries if they 
fulfill certain conditions, in particular stationarity and ergodicity. The method was 
applied to diffusion with and without nonlinear microscopic boundary conditions and 
to porous media flow. Some examples for stochastic geometries where given for solid 
and porous microstructures. Note that we did not treat homogenization of problems 
dfU^ — div (e^D^Vw^) = / with e^-scaled diffusion. However, such problems can be 
treated even easier than the example calculation above with results of the form dtUo — 
div^ (DVujUq) = f. 

The results of this article can be considered as one more step towards understanding 
homogenization in non-periodic heterogeneous media. A big advantage of the presented 
approach is, that formal calculations require only few new theory. The only formal 
difference in the calculation is, that Y has to be changed into fl and y into u. Although 
the mathematical theory behind the rigorous calculations gets much more complex, 
the method's user is not bothered by that. Indeed, besides some critical points in the 
calculations which were pointed out, one may not even think of the fact that one is 
doing stochastic calculations. The limit equations are independent on the particular 
realization which was used for homogenization. This is a feature of ergodicity and 
stationarity of the geometries. 

A major issue for further investigations is the search for suitable models for geometries 
of natural heterogeneous media. Also, from the mathematical point of view, it would 
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be nice to develop a "stochastic unfolding" corresponding to the periodic unfolding as a 
further theoretical foundation of the presented asymptotic expansion. Finally remark, 
that the intention of this paper was not to give a rigorous introduction to stochastic 
homogenization or the mathematics behind it but only to demonstrate that for non- 
mathematicians, switching from periodicity to stationary ergodic geometries can be 
easily achieved and would not bother the calculations or the results. 



Appendix A. Stochastic Geometries 

This section follows [8j to introduce some basic mathematical concepts on random 
geometries. In particular, random closed sets and random measures will be defined 
rigorously. The connection to the periodic case was explained in detail in [8] and is 
based on similar ideas as the connections between periodic and stochastic asymptotic 
expansion. It will be shown that the assumptions on A, B and F in section 3.1 are 
reasonable and can be made w.l.o.g.. For more information on random closed sets, 
refer to Matheron |14| or Molchanov |17| . 

Let J-'(]R") denote the set of all closed sets in M". Then for K G compact and 
V gW^ the following sets can be defined: 

(A.l) J^v ■■= {F G J^(M") I FnV y C M" open set 

(A.2) := {F e J^(M") I FnK = K compact set 

The Fell topology on J-'(M"') is created by the sets J-'y, J-"^ for all open V and compact K 
and according to [T^, (J-'(M"'), Tj-) is compact, Hausdorff and separable. The Matheron- 
a-field ajr is the Borel-cr-algebra created by the Fell-topology. 

For a probability space {Q,a,fi), a Random Closed Set (RACS) is a measurable map- 
ping 

A : (fi, a, /i) — > (J", ajr) . 
It is the aim of this short introduction to demonstrate that such RACS have the nice 



properties which we assumed on A, B and F from section |3.1[ Remark that up to now, 
a RACS is a mapping from a probability space into the set of closed subsets of M". 

In what follows, M.{W^) respectively M. denotes the set of all locally finite Borel 
measures on M". In Particular, C denotes the Lebesgue measure and the m- 
dimensional Hausdorff measure. The smallest topology such that 

M^R, /i I — > J fdjl 

is continuous for all / G Co(M") is in general called the Vague-topology 7v on A^(]R"). 
The Borel-(T-field of this topology is denoted by Uv respectively by B{M.). 

Theorem 1. [6j (A^,7v) is a separable metric space. The a-algebra cry on Ai created 
by the Vague topology 7v is the smallest one, such that fj, i-^ /^(-B) is measurable for 
every bounded and measurable set B C M". 

A random measure on M" is a measurable mapping {Q,a,fi) — )■ {A4,a^), u i— )■ fi^^, 
where {Q, a, /x) is a probability space. Random closed sets and the random measures 
are related due to the following Lemma, which implies that a random closed set always 
induces a random measure. 



Lemma 2. ([21] Theorem 2.1.3 resp. Corollary 2.1.5) 
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Let Tm C J-" be the space of closed m-dimensional sub manifolds of such that the 
corresponding Hausdorff measure is locally finite. Then, the a-algebra ajr fi J-^ is the 
smallest such that 

M]^ : j-^ M r W^ij: n b) 

is measurable for every measurable and bounded i? C M". 



It was part of the argumentation in [8] that due to Lemma |2] for any initial random 
closed manifold A : cr, /x) — )■ Tm the assumptions ^ d M. and a = B{J^) fl Q can 



be made w.l.o.g.. We introduced dynamical systems in section 3.1 and also defined 
ergodicity. However, one normally assumes only measurability of r but the continuity 
is a rather direct consequence of C ^Sj. 

A random measure is said to be stationary if for every u & Q: T{x)fiuj = P-r^ui with 

T(xK(5) :=/i,(5 + a;) V5 G 

and a RAGS A is called stationary if XA{Tyui){x) = XAiu)){x + y) where Xa{ui) is the 
characteristic function of A{u). This is slightly different introduction of stationarity 



than in section 3.1 However, we will see below that these definitions of stationarity and 



ergodicity already guaranty that RAGS have the properties which we claimed in section 



3.1 First, it is necessary to show that there is an equivalent of a Hausdorff-measure 
on r. This Hausdorff measure will be fir,p, stated by the following 

Theorem 3. (Mecke [151 16]: Existence of Palm measure) 

Let C be the Lebesgue-measure on with dx := dC{x) and {Q,a,fi) as above. Then 
there exists a unique measure fi-p on Q such that 



(A.3) / / f{x,Txu)dfiui{x)diJ,{u) = / / f{x,u)dfiv{u)dx 

Jn Jr" Jn 

for all i3(]R") x i3(fi)-measurable non negative functions and all /ip x C- integrable 

functions. 

fi-p as in Theorem [3] is called the Palm measure of /i^. Note that ^ is the Palm measure 
of C For a random closed manifold F : {Q,a,fi) — )■ J-^ u h-)- F(cj) the random 
Hausdorff measure will be denoted hj u ^ f^r{uj) and the Palm measure by fir,v- 

Theorem 4. (Ergodic Theorem [6j) 

Let the dynamical System be ergodic and assume that the stationary random mea- 
sure iJ,^] has finite intensity. Then 

(A.4) lim } I g{T^u)dfiuj{x) = [ g{u)dfi-p{u) 



t—^oo 



t"£(A) 



tA Jn 



for fl almost surely in u for all bounded Borel sets A D {|x| < 1} and all g G L^{Q, fi-p) 

For /ip(.^^(i?) := e"'fir{u}){s^^B) the latter theorem yields after a transformation of 
variables 

(A.5) lim / f{x,Tiiuj)dfi'l,,.= / / f{x,uj)dfir,vdx. 

^'^^Jq ' JqJn 

for / = (I){x)iIj{uj) with (p G C^{Q) and G Cb{^) and therefore also for / G 
C^(Q; Cb(fi)). This is a first indicator why we may consider fir,r as a Hausdorff 
measure. It will be more striking by the following 

Lemma 5. [8j There is a measurable set, also denoted by F C fl, with Xr{uj){x) = 
Xt{tx(^) for C + /ir{(^)-almost every x for /^-almost every u. Furthermore /i(F) = and 
lirA^\T) = 0. 
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The same proof also provides such a characteristic function for a n-dimensional random 
closed set A. F := dA is also a random closed set, which can be verified using J-y from 



(A.l) with V = A\dA. In case T is regular enough, there are subsets A,T C 
such that XA{Lu)i^) ~ Xa{'Tx<^) and Xr{uj)i^) ~ Xr{'Tx<^)- This is equally possible for 

B := ^M^y^^ U f . Thus, the assumptions made on A, B, and T made in section 



3.1 are now justified. Note that it is also possible to define i/p : F — )■ M" such that 



^r{uj){x) = h'ri'Tx'^) and ur thus can be considered as normal field on F. 



Note that in section 3.1, we started with A,B,T C and defined A(uj), B{lo) and 
to avoid these rather technical preliminaries. We conclude the mathematical part 
by some remarks on function spaces on as well as on Gauss' theorem and Poincare 
inequalities. 



Continuity and differentiability of functions were introduced in section 3.1 It is of 
course also possible to define arbitrarily high differentiability C^{Vt). Since f2 is a 
separable metric space equipped with a measure yU, it is possible to define 

ViVt^ cr, /i) := |/ : — 7- M : / is o" — measurable, j \ f\^ (i/i < ooj- 

These L^-spaces have good properties, in particular they are separable and for any 
finite measure /i on {VL,B{VL)) exists a countable dense set of C^(f2)-functions in 
LP(fi,i3(fi),/i) (1 <p < oo, /i cT-finite )[8]. 

It is of course possible to define various kinds of Sobolev spaces on Vt. The reader is 
referred to |2Sl [271 E] • 

Finally, in order to justify some of the calculations in section [4], it is vital to proof a 
Gauss-like theorem on stochastic spaces as well as some Poincare inequalities. 

Lemma 6. iSj For all G CI{Vl) and <p G Q^fi)" holds: 

(A.6) / V'Vf^d/i = - / {V^)(j)dfi 

Jn Jn 

Since the technique is very often used, we give the abbreviated proof from [8]: 

Define Qm ■= [—m,m]"'. Then, since (p and are bounded, the Ergodic Theorem |4] 

leads to: 



Jn '"^"^ (2m)" Jq^ 

1 

(2m)" 



lim \ ( - / {Vtp){T^u)(f){T^u)dx 



m—^oo 



iyQjx){<p^)iT,u)d'H''-\x: 

dQm 



(2m)^ 
{V^)(t)d^i 



Qn 



n 



where the integral over the boundary vanishes due to the fact that the surface of 
grows with (2m)"~^.QED 

Note that ( |A.6 ) would also hold for an integral over A if -i/; = on F or V^/^ ■ z/r = on 
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From the corresponding Poincare inequalities on M", one may conclude in the same 
way 

[ + < c( [^ + [ ivijf) y^eclin) 

Jn \Jn Jn J 

[ {ip^ + \Vipf) < C I |V^|^ ^i) E Cl{A) with iP = OonT 

J A J A 

where we have to assume that the constant C for the realizations is independent on u: 

jQmr\A(uj) JQmr\A{u)) 

Once such inequalities are established for C^-functions, it is easy to expand them to 
corresponding Sobolev spaces on Vl. (See P, [Hj for more complicated examples) 
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